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I. INTRODUCTION 



The principle of factorization underlies all theoretical predictions for hadronic processes. 
Simply put, factorization is the statement that short and long distances contributions to 
physical processes can be separated, up to corrections suppressed by powers of the relevant 
large scale in the process. The predictive power gained from this result stems from the fact 
that the incalculable long distance effects are universal, defined in an unambiguous way in 
terms of matrix elements. As a consequence, the non-perturbative long distance effects can 
be extracted in one process and then used in another. In general, proving factorization is a 
difficult task |lj. The proof of factorization in Drell Yan processes, for instance, took several 
years to sort out [[| (for reviews on factorization see [J). Indeed, there are still some 
processes such as B — > -kit where a proof of factorization only exists at one-loop ||. 

Given that we would like to retain our predictive power over the largest possible range of 
energies, we are compelled to understand power corrections to the factorized rates. These 
corrections are not necessarily universal, and as such, the relevant size of the power cor- 
rections are process dependent. In processes for which there exists an operator product 
expansion (OPE), there is a systematic way in which to include power corrections. How- 
ever, for a majority of observables we do not have an OPE at our disposal, and the nature 
of the power correction is not always known. For instance, in the case of shape variables 
there is still some on going discussion about the form of subleading corrections |7], ||. 

The purpose of this paper is to show that an effective theory framework can be used 
to simplify proofs of factorization and describe processes with an operator formalism. To 
do this we extend the soft-collinear effective theory (SCET) developed in Refs. || |H], [11 



12] , |I3[ |, to high energy processes. It should be emphasized that there are several other 
useful advantages in using an effective field theory (EFT). For instance, the EFT makes any 
symmetries which emerge in the Q — > oo limit manifest in the Lagrangian and operators, 
and allow statements to be made to all orders in perturbation theory. The calculation 
of hard coefficients reduces to simple matching calculations, where subtracting the EFT 
graphs automatically removes all infrared divergences from the QCD calculation. Perhaps 
most importantly, it provides a framework for systematically investigating power corrections. 
Finally, the EFT framework allows standard renormalization group techniques to be used 
for the resummation of logarithms that are often necessary in calculating rates for certain 



high energy scattering events || [10], [14], [15[]. The factorization formulae that we prove in 
this paper are not new, but serve to illustrate our approach in familiar settings. The results 
are valid to all orders in a s and leading order in the power expansion. The simplicity of our 
approach lies in the fact that factorization occurs at the level of the SCET Lagrangian and 
operators, and is facilitated by gauge symmetry in the EFT. This provides the advantage 
that our proofs do not rely on making use of Ward identities and induction, or on specifying 
a particular gauge. 1 Furthermore, it becomes rather simple to derive factorization formulae 
for a myriad of processes, since many results are universal. The examples given here serve 
to illustrate these simplifications. Developments on the issues of power corrections and 
resummations are left to future publications. 2 

In section 11 we review the construction of the SCET. The formalism developed in Refs. || 



1 In fact our factorization proofs rely heavily on the gauge symmetry structure of SCET. When a gauge 

fixing term is required for explicit calculations we use general covariant gauges. 

2 For recent work on subleading corrections in SCET for heavy-to-light transitions see Ref . Co] . 
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TO] , |TT| , |l"2"f is extended to include two types of collinear particles moving in opposite directions 
in section LIB , and factorization for 7* to two collinear states is discussed as an example. In 
section [11 C we define the non-perturbative matrix elements such as the parton distribution 
functions that will be needed for the processes presented in the paper, and in section |11 D 



we 



discuss some of the symmetries in SCET that may be used to place restrictions on matrix 
elements. In the remaining sections we give various examples on how factorization theorems 
emerge in the effective theory language. In section |TTT] we prove factorization theorems for 
two exclusive processes, namely the n-j form factor, and meson form factors (7*M — ► M) 
for arbitrary spin and isospin structure. In section [TV] two inclusive processes are treated, 
namely DIS (e~p — > e~X) and Drell-Yan (pp —>■ X£ + £~), and we also give results for deeply 
virtual Compton scattering (7*^ — > 7^p). In these processes we include all leading power 
contributions in the factorization proofs (even if the operators are only matched onto at 
higher orders in perturbation theory such as for the gluon distribution functions). Our 
conclusions are given in section [V]. In appendix |A] we show how auxiliary fields can be used 
to prove the simultaneous factorization of soft fields from collinear fields for particles in 
back-to-back directions. 



II. FORMALISM 

Effective field theories provide a simple and elegant way of organizing physics in processes 
containing widely disparate energy scales. In constructing an EFT, some degrees of freedom 
are eliminated, and the remaining degrees of freedom must reproduce all the infrared physics 
of the full theory in the domain where the EFT is valid. The EFT is organized by an 
expansion in A, defined as the ratio of small to large energy scales. As a useful guideline 
the following steps are used to identify the infrared degrees of freedom: 1) Determine the 
relevant scales in a problem from the size of the momenta and masses of all particles that 
can make up the initial and final states, 2) Construct all momenta from these scales whose 
components correspond to propagating degrees of freedom, and which have offshellness less 
than the large scale, i.e. p 2 — m 2 < Q 2 . Effective theory fields are then constructed for each 
unique set of these momenta. 

We will be interested in an EFT with particles of energy Q much greater than their mass. 
The dynamics of these particles can be described by constructing a soft-collinear effective 
theory (SCET). This theory is organized as an expansion in powers of A ~ P±/Q, and offshell 
fluctuations with p 2 ^> (QX) 2 are integrated out. In section |11 A| we begin by describing this 



procedure and comparing the construction to other EFT's. We then give a brief review of 
the soft-collinear effective theory developed in Refs. || [10, [II], [T2|, [H| . We do not attempt 



to give a comprehensive treatment, but instead emphasize the main results and refer the 
reader to the literature for details. In section |II B| we extend the formulation of SCET to 
describe processes with collinear particles moving in back to back directions, and prove the 



factorization formula for 7* to two collinear states as an example. In section [II C| we define 
the non-perturbative matrix elements that are needed for our examples, then in section |IID 
we discuss some of the symmetry properties of collinear fields and currents. 
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A. Soft-Collinear Effective Theory 



In the standard construction of an EFT one removes the short distance scales and massive 
fields by integrating them out one at a time. A classical example is integrating out the W 
boson to obtain the effective electroweak Hamiltonian with 4-fermion operators. However, 
in some situations we are interested in integrating out large momentum fluctuations without 
fully removing the corresponding field. The simplest example of this is Heavy Quark Effective 
Theory (HQET) fPT| , which is constructed to describe the low energy properties of mesons 



with a heavy quark. Here the heavy anti-quarks are integrated out and only heavy quarks 
with fluctuations close to their mass-shell are retained. This is accomplished by removing 
fluctuations of order the heavy quark mass rriQ with a field redefinition |Hf 



ij(x)=J2e- imQV - x h v (x), (1) 

V 

where v is the heavy quark velocity and h v is the field in the EFT. While d^ip{x) ~ itiq i/j(x), 
the effective field has d^h v (x) ~ Aqcd^(^), indicating that it no longer describes short- 
distance fluctuations about the perturbative scale ttiq. Instead these effects are encoded in 
calculable Wilson coefficients. The HQET degrees of freedom with offshellness p 2 ~ Aq CD 
are the heavy quarks, soft gluons, and soft quarks. 

Similarly, for collinear particles with energy Q ^> m, one needs to remove momentum 
fluctuations ~ Q while retaining effective theory fields to describe smaller momenta. How- 
ever, unlike heavy quarks the collinear particles have two low energy scales. Consider the 
light-cone momenta, p + = n ■ p and p~ = n ■ p where n 2 = n 2 = and n ■ n = 2. Here n 
parameterizes a light-cone direction close to that of the collinear particle and n the opposite 
direction (eg. for motion in the z direction = (1,0,0,1) and = (1,0,0,-1)). For a 
particle of mass m < p± <C Q, we have p~ ~ Q, and a small parameter A ~ P±/Q- The 
scaling of the p + component is then fixed by the equations of motion p + p~ + p\ = m 2 , so 
that (p + ,p~,p±) ~Q(A 2 ,1,A). 

The appearance of two small scales, QX 2 <C QX Q, is similar to the situation in 
non-relativistic QCD (NRQCD), which is an EFT for systems of two heavy quarks with 
an expansion in their relative velocity j3. In a non-relativistic bound state the momentum 
of a heavy quark is p ~ ^qP, but the equations of motion E = p 2 /(2mg) make the 
energy E ~ toq/? 2 , giving scales rriQf3 2 -C mQj3 -C rriQ. The two low energy scales can be 
distinguished by following Eq. (0) with a further field redefinition [|Tj| h v (x) = J2 P e ip ' x V ; p(^), 
so that derivatives on ip p only pick out the m(3 2 scale. The on-shell degrees of freedom are 
then the heavy quarks, soft quarks and gluons with p 2 ~ (rriQv) 2 , and ultrasoft quarks and 
gluons with p 2 ~ {rriQV 2 ) 2 . 
SCET fields: 



For collinear particles the analogous field redefinitions are |H1 ITT 



^ = E E e " ip ' x< M x ) ' (2) 

n p 

where the collinear fields <p n tP are labelled by light-cone vectors n and label momentum p. 
Here p contains the n ■ p ~ Q and p± ~ QX momenta so that <9 M 0„ iP ~ (QA 2 )0 njP . The 
field 4> p can either be a quark or gluon field. Similar to if) p , the missing ~ Q fluctuations 
are described by Wilson coefficients and the ~ QX labels simplify the power counting by 
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distinguishing the QX and QX 2 scales. Now 



so collinear particles and antiparticles are contained in the same effective theory field, but 
have momentum labels with the opposite sign. In the large energy limit the four component 
fermion spinors contain two large and two small components. One therefore defines collinear 
quark fields £ niP which only retain the large components for motion in the n direction and 
satisfies ^£ n ,p = 0. For these fields Ct, P /Cn, P destroy /create the particles/antiparticles with 
large momentum n ■ p > |TT[. For collinear gluons = A^_ q , and (A£^ q ) + / (A£ q )~ 



destroy /create gluons with n-q > 0. 

For simplicity we will ignore quark masses and only consider massless u and d quarks. 
For the processes considered here SCET then requires three types of degrees of freedom: 
collinear, soft, and ultrasoft (usoft) fields. These are distinguished by the scaling of the 
light cone components (p + ,p~,p ) of their momenta: (A 2 ,1,A) for collinear modes in the 
n direction (A nq , £„ )P ), (A, A, A) for the soft modes (A q , q^), and (A 2 , A 2 , A 2 ) for the usoft 
modes (A us , q us ). The soft modes are labelled by their order QX momenta, so A s q and q^ 
are essentially just momentum space fields. The usoft fields have no labels and depend only 
on the coordinate x. The fields are assigned a scaling with A to make the action for their 
kinetic terms order A || [12). For instance ~ A, A^ q ~ (A 2 ,1,A), A q ~ A, and 
Aus ~ At leading order only order A vertices are necessary to correctly account for all 
order A Feynman diagrams. 

In HQET only external currents with momenta of order rrib can change the label v. Thus 
the Lagrangian has a superselection rule forbidding changes in the four- velocity of the heavy 
quark [|17[ 18f . In NRQCD the v labels are also conserved, but the smaller momentum labels 



p are changed by operators in the effective theory such as the Coulomb potential. A novel 
feature of SCET is that interactions in the leading action can change both the large and 
small parts of the momentum labels p^. However, only external currents can change the 
direction n of a collinear particle, so this label is conserved. Thus, for each distinct direction 
n a separate set of collinear fields are needed. In the remainder of this section we will 
restrict ourselves to collinear particles with a single n. We will generalize the discussion to 
the case of two back-to-back directions and discuss the factorization of collinear particles 
with different n's in section I11B. 



Since in SCET interactions can change the order Q label momenta it turns out to be 
very useful to introduce a label operator, V 1 [|TT|, for which the collinear fields satisfy 
'P^£,n,p = More generally, acts on a product of labelled fields as 

f(vn U ■ ■ ■ o Pi • • • J = f(pi+. u . . . 0pi . . . j . (4 ) 



qi rq m ^Pl Vp n J J \F1 i • • • i Fn HI • • ■ Hm) \Tq\ Yq m rVl 

so conjugate field labels come with a minus sign. The operator acts to the right, while 
the conjugate operator acts to the left. As explained in Ref. the label operator allows 



all large phases to be moved to the front of operators with a factor exp(— ix-V). This phase 
and the label sums can then be suppressed if we impose that interactions conserve label 
momenta and that the momentum indices on fields are implicitly summed over. Basically, 
for labels p and p' and residual momenta k and k' 

J d 4 X e i(P'-P+k'-k)- X = 5{p _ p /) J e i{k'-k). X ? (5) 
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so that the label and residual momenta are individually conserved. (Although technically 
the label momenta are discrete we abuse notation and use 5(p—p') rather than 5 PtP i because 
it makes the subscripts easier to read.) For convenience we define the operator V to pick 
out only the order A labels on collinear fields, and the operator V 1 to pick out only the 
order A labels. For the matrix element of any collinear operator O, momentum conservation 
constrains the sum of field labels ]TT[] , giving 



(M ntPl \[f(V)0]\M n , P2 ) = f(n.(p 2 - Pl )) (M n , Pl \0\M n , P2 ) , (6) 
for any function /. 

For a single n the Lagrangian can be broken up into three sectors: collinear, usoft, and 
soft. We therefore write 

where we have made the field content of each sector explicit. We will discuss each of these 
terms separately. 
Collinear sector: 

As explained in detail in Ref. gauge invariance in SCET restricts the Lagrangian 
and allowed form of operators. Only local gauge transformations whose action is closed on 
the effective theory fields need to be considered. These include collinear, soft, and usoft 
transformations. Each of these vary over different distance scales, with collinear gauge 
transformations satisfying d^U n (x) ~ Q(X 2 , 1, A) U n (x), soft satisfying d^V s {x) ~ QX V s (x), 
and usoft transformations with d^V us (x) ~ QX 2 V us (x). All particles transform under V us (x) 
and usoft gluons act like background fields for collinear particles. Invariance under U n (x) 
requires a collinear Wilson line built out of the order A gluon fields [TO, [TTj 



WJx) 



E exp( -g -=n-A nq {x)) 
perms V y J 



(8) 



Here the operator V acts only inside the square brackets, the n on W n refers to the direction 
of the collinear quanta, and W n is local with respect to x (corresponding to the residual 
momenta). Taking the Fourier transform of 5(u — V)W n (0) with respect to uo gives the 
more familiar path-ordered Wilson line W n (y, — oo) = P exp[ig ds n-A n (sn)\ . Under a 
collinear gauge transformation W n transforms as W n — > U n W n . An invariant under collinear 
gauge transformations can therefore be formed by combining a collinear fermion and 
the Wilson line W n in the form 



Wl(x) £ n>p (x) . (9) 

This combination still transforms under an usoft gauge transformation, W^ ntP — > 
V us (x) W^ ntP . We will often suppress the x dependence of the combination W\ 

Integrating out hard fluctuations gives Wilson coefficients in the effective theory that are 
functions of the large n-pi collinear momenta, C(n-pi). However, collinear gauge invariance 
restricts these coefficients to only depend on the linear combination of momenta picked 
out by the order A operator V [|T_TJ. In general the Wilson coefficients are then functions 
C(V,V^) which must be inserted between gauge invariant products of collinear fields. In 
general the Wilson coefficients also depend on the large momentum scales in a process such 
as Q and the renormalization scale /i. 
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To construct the collinear Lagrangian one can match full QCD onto operators with 
collinear fields that are invariant under usoft and collinear gauge transformations. The 



collinear Lagrangian at order A is [TO] , [TT] , |T2 
r 



1 



2g 2 



tr 



iV» + gA% jq ,iW + gAl al 



n,q' 



(10) 



where £f are gauge fixing terms, i_D M = i<9 M + gA^ s , and 



nr - 



= — V + Vt+'—in-D 



(11) 



Since usoft gluons act as background fields in the collinear gauge transformation the cou- 
plings, g(yu), for both types of gluons must be identical. 
Usoft and Soft sectors: 

The usoft and soft Lagrangians for gluons and massless quarks are the same as those 
in QCD. From Eq. (^) we see that collinear quarks and gluons interact with usoft gluons, 
however at order A only the n-A us component appears in Eq. (pi]). In order to prove 
factorization formulae it is essential to disentangle the collinear and usoft modes. This can 
be done by introducing an usoft Wilson line 



YJx) 



P exp (^ig J 



ds n-A^Jsn 



(12) 



where the subscript n on Y n labels the direction of the Wilson line (we emphasize that this is 
different from the meaning of the subscript on W n in Eq. (§)). An usoft gauge transformation 
takes Y n -> V us Y n . In Ref. @] it was shown that the field redefinitions 
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i hp 



1 n Sn,p 



V • 



A 1 * 

n,p 



Y 4(°)/ i vt 
1 n -^ritP 1 n i 



(13) 



imply W n = Y n Wn^Y^ and decouple the usoft gluons from the collinear particles in the 
leading order Lagrangian 



£c,n [£n,p) Af^ „, ?1."v4 us ] 



r r^(°) 4(°)a» nl 



(14) 



Thus, the new collinear fields with superscript (0) no longer interact with usoft gluons or 
transform under an usoft gauge transformation. Since the field redefinitions do not change 
physical S matrix elements, the new fields give an equally valid parameterization of the 
collinear modes. The leading SCET Lagrangian therefore factors into separate collinear 
and usoft sectors. This alone does not guarantee factorization in operators and currents, 
since after the field redefinition these operators may still contain both usoft and collinear 
fields. However, the field redefinition makes factorization transparent since identities such as 
Y^Y n = 1 may be applied directly to the operators. This will become clear in the examples 
in sections |TTTj and |TV| . 

The coupling of soft gluons to collinear particles differs from the usoft-collinear inter- 
actions. Interactions of a soft gluon with a collinear particle results in a particle with 
momentum p ~ Q(A, 1, A), so soft gluons can not appear in the collinear Lagrangian. These 
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offshell particles have p 2 ~ Q 2 X and since Q 2 X 3> (QX) 2 these offshell quarks and gluons 



can be integrated out. At leading order in A it was shown in Ref. [12] that in operators with 
collinear fields this simply builds up factors of a soft Wilson line S n involving the n ■ A s 
component of the soft gluon field, 



Sn 



exp 

perms 



n ' As ' q 



(15) 



The factors of S n appear outside gauge invariant products of collinear fields, and their 
location is restricted by soft gauge invariance. 



B. SCET for n and n collinear fields 

In this section we extend SCET to include the possibility of collinear fields moving in 
different light-cone directions: n 1; n 2 , n 3 , . . . . These directions can be considered to be 
distinct provided that rii-rij ^> A 2 for i =^ j. This follows from the fact that if n\ • rig ~ A 2 
then the directions n\ and n 2 are too close to be distinguished. For example, a momentum 
P2 = Qn 2 can be considered to be collinear in the n\ direction if n\-p 2 = Qn\-n 2 ~ QX 2 , since 
this is the correct scaling for the small momentum component of an ni-collinear particle. 

For simplicity we will only consider the case of back-to-back jets corresponding to collinear 
particles moving in the n and n directions. These are clearly distinct since n-n = 2. Collinear 
particles in the n direction have (+, — , _L) momenta ~ Q(l, A 2 , A), and the n - p ~ 1 and 
p± ~ A momenta appear as labels on the corresponding fields: £ s>p and A£ p . Emission of 
a collinear particle moving in the n direction from a collinear particle in the n direction 
results in a particle with momentum k ~ Q(l, 1, A) and offshellness k 2 ~ Q 2 . These offshell 
modes are integrated out to construct the SCET, so collinear modes in the n direction 
do not directly couple to collinear modes in the n direction. A distinct set of collinear 
gauge transformations is associated with each of n and n, and fields in one direction do not 
transform under the gauge symmetry associated with the opposite direction. Two Wilson 
lines W n (x) and Wn(x) are necessary (defined as in Eq. (§)), and they appear in a way that 
makes collinear operators gauge invariant. For instance the combinations 

W^p, Wk n , p (16) 

are invariant under collinear gauge transformations in the n and n directions, respectively. 
We also require two types of label operators, V as before, and an operator V to pick out n-p 
labels that are order A . Thus, V and V act only on the n and n collinear fields respectively. 
(The label operator V^ 1 still picks out order A momentum components and therefore acts 
on both n and n fields.) With two collinear directions, decoupling usoft gluons requires 
introducing both Y n and Y n Wilson lines, defined as in Eq. (|P2D, but along the n or n paths 
respectively. Finally, integrating out ~ Q 2 X fluctuations at leading order induces both S n 
and S n soft Wilson lines defined analogous to Eq. (|TJ). This is discussed in greater detail 
in Appendix |A| where we show explicitly to all orders in g that integrating out the Q 2 X 
fluctuations causes 

WlCn, P ^ S n W^ n>p , Cn, P W n - ln, P W n Si , 



8 



Relations for operators with collinear gluon fields are also derived in Appendix [A[ 

Note that we have not included "Glauber gluons" with momenta p M ~ (A 2 , A 2 , A), which 
are kinematically allowed in t-channel Coulomb exchange between n and ft collinear quarks. 
In determining the relevant degrees of freedom we have assumed that Glauber gluons are 
not necessary to describe the infrared for the processes considered in this paper. Intuitively, 
this can be seen from the fact these gluons are instantaneous in both time and longitudinal 
separation, and only could contribute when the n and ft jets overlap for a duration of order 
1/(QA 2 ) in a space-time diagram. In processes with a hard interaction the overlap scale 
is always much shorter than this (however this need not be the case in processes such as 
forward scattering). For the Drell-Yan process more quantitative arguments can be found 
in Refs. |§. 



At order A it is not possible to construct a gauge invariant kinetic Lagrangian with 
terms that involve both n and ft fields. Thus, the n and ft collinear modes are described 
by independent Lagrangians (however n and ft modes may still both appear in an external 
operator). The collinear sector of the SCET Lagrangian is therefore 

-£-c,n[£n,p) ■^■n,gi n ' A U s] 4~ Cc,n [£ra,p; -^-n,qi n ' A us ] . (18) 

Making the field redefinitions 

t y t(o) Afj, _ y a{o)hy\ 

Sn,p 1 n S,n,p > n,p 1 ri^n^p 1 n j 

c _ y tip) a» _ y_^(_°)^y_t (ig) 

S>n,p 1 n Sn,p > ■ rL n,p 1 n r >-ri,p 1 n j \ ±,J J 

gives W n = YnWn^Y^ W n = Y n W^Y^, and usoft degrees of freedom once again decouple 
from the collinear modes since 

£c,n [£,n,pi -^n,qi n ' A U s] 4" £c,n [£,n,pi ^-n,qi n ' A U s] £c,n [£,n)p , > ^n,g^' 0] 4"/^c,n[^n,p) ^-n,q i 0] • (20) 

Thus, usoft gluons are removed from the collinear Lagrangian at the expense of inducing Y n 
and Y n factors in operators with collinear fields. In certain cases the identities Y^Y n = 1 and 
Y^Y n = 1 can be used in these operators to cancel usoft gluon interactions. Perturbatively 
these cancellations would occur by adding an infinite set of Feynman diagrams. 

To see in more detail how this works consider the simple example of the 7*-production 
of back-to-back collinear states X n and X n . The full theory current $(x)Y r ip{x) matches 
onto an effective theory operator O nn - Naively one might guess that the SCET operator 
mediating this process is 

O n n ^n,piF^fi,p2 ■ (21) 

However, this operator is not invariant under the collinear gauge transformations U n and 
U n , so the process is instead mediated by the invariant operator 

O n n = Up,W n TWUn,p 2 . (22) 

A hard matching coefficient C(V, V\ V, V') can be inserted in any location in the operator 
that does not break apart the gauge invariant combinations of fields in Eq. ([16]). The 
operators V and V in the coefficient only pick out momenta that are order A in the power 
counting. Thus, V does not act on fields in the n direction and V does not act on fields in 
the ft direction, and the most general result is 

Onn = U Pl W n T C{V\ V)Wi in,p 2 • (23) 
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Next, we integrate out the offshell Q 2 \ fluctuations which induces additional soft Wilson 
lines in the operator. This is discussed in detail in Appendix [X| and from Eq. (|T?D gives 

o nn = L, P ,w n si r c{v\ v)s n wt £ n>P2 . (24) 

Note that V and V do not act on the fields in the soft Wilson lines since soft gluons carry 
only order A momenta. Finally, we can make the usoft gluon couplings explicit by switching 
to the (0) fields using Eq. ([19]): 

Onn = &lM° )Y nSi r C(P\ V) S n YM 0) ^n°l 2 • (25) 

This operator is manifestly invariant under collinear gauge transformation in the n and n 
direction, as well as under soft and usoft gauge transformations. 

To separate the short distance Wilson coefficient from the long-distance operator one 
introduces convolution variables u and u' to give 

O nn = Jdujduj'C(u;,u')O n n(uj,uj') , (26) 

o nfl (<"V) = fe^ ^^ 1 -^)^^^^^^-^)^^ ^^ 



The function C(u, u/) contains all the short distance physics and is determined by matching 
the full theory onto this effective theory operator. O nn {uj,uj') contains all the infrared long- 
distance QCD contributions at leading order in A. 

Now consider the matrix element of the production current between (X n X n \ and the 
vacuum. Taking the 7* to have large time-like momentum g M = (Q, 0, 0, 0) (and zero residual 
momentum) we have 



d 4 Xe~^ X (X n X n \tfj(x)T^(x)\0) = jd 4 x(X n Xn\Onn(x)\0) 

d A xe ikx ((X n X n )(k)\Onn(x = 0)|0) = <(X n X s ) (0) \O m (0) |0> . (27) 



In the first step the conservation of the large label momentum q was made implicit in 
the matrix element (c.f. Eq. (El)). Since we are in the center-of-mass frame the X n has 
large momentum n ■ p = Q, and the X n has momentum n ■ p' = Q. Now using translation 
invariance, we see that the remaining x integral forces the \X n X n ) state to have zero residual 
momentum. Using Eq. (EBT) this matrix element is equal to 



duj duj'C{uj, cu / )(X n X fi |O n f i (a;, a/) |0) = jduj doj' C(lu,lu') J n (co>) T S nn Jn(^') , (28) 



where the functions C, J n , J n , and S also depend on the renormalization point fi. Here we 
have used the fact that both X n and X n can be described entirely by collinear particles in 
the n and n directions respectively. Since the Lagrangians for the collinear, soft, and usoft 
fields are factorized the remaining matrix element splits into distinct matrix elements for 
each class of modes. These matrix elements are 

J n {u) = {XjQ/2)\^W®8(V*-u>)\Q), 

Uu') = (X n (Q/2)\S(V-u;')W^knt\0) , 

Snn = (0\Y^SlS n Y n \0) , (29) 
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(and are matrices whose color, spin, and flavor indices are suppressed). Note that J n , J n , and 



S n fi are explicitly invariant under the collinear, soft, and usoft gauge transformations [12 
of SCET, but still transform globally under a color rotation. Now using the momentum 
conservation identity in Eq. (|6|), the large momentum of the X n and X n states set V — > Q 
and V> — > Q. Label conservation also implies that the total perpendicular momentum of 
each of J n , J n , and S nn is zero. The sum over u and uj' can then be performed to give the 
final factorized form 

C(Q,Q) J n (Q)TS nn J n {Q). (30) 

Although rather idealized, the above example illustrates the main steps needed to derive a 
factorization formula. Taking X n and X n to be single quark states the result in Eq. (]30f) also 
agrees with the factorization formula for the onshell production form factor for qq [[H], |22|[. 3 
In the above example the factors of S and Y in the operator in Eq. fl25|) do not cancel. In 



the examples we will consider in sections III and [V] there are several operators at leading 



order, however the factors of S and Y cancel in observable matrix elements. The collinear 
matrix elements of long-distance operators, such as those in Eq. (|29"D, are the ones that have 
interpretations as structure functions or wave functions. In the next section we give the 
operator definitions for these functions that will be needed in the remainder of the paper. 



C. Non-perturbative Matrix Elements 

Predictions for hadronic processes depend on universal matrix elements that are not 
computable in perturbation theory. For exclusive processes these include light-cone wave- 
functions and form factors, while for inclusive processes they include parton distribution 
functions and fragmentation functions. In this section we define matrix elements in SCET 
that are needed for our examples. All the collinear operators considered here decouple from 
usoft gluons since they are local in the residual coordinate x and because Y\x)Y(x) = 1. 
Thus 

U P >WTW^ n , p = ^l> wi ° )TWm Znl , (31) 

and expressions with and without the (0) superscript are equal. For convenience we will 
write the collinear fields without the superscript in the remainder of this section. 

Consider first the light-cone wavefunctions. For the pion isotriplet ir a , the wavefunction 
4> n (x) is conventionally defined by ^ZE 



(7r a (p)\^(y)Y^^=Y(y,x)^(x) |0> = -ifJ a Y [ dz e^^ 1 '^ ^(/i, z), (32) 



o 



3 In this case depending on the choice of infrared regulator(s), it may not be possible to distinguish the Y n 
and S n Wilson lines in S n n{^)- For instance if one chooses Air ~ QX then the usoft gluons give scaleless 
loop integrals and can be dropped, so that S^SnYn — > S^Sn. If instead one chooses Air ~ QX 2 then 
the soft gluons give scaleless loop integrals (they simply act to pull-up the ultraviolet divergences in the 
usoft integrals to the hard scale [^3[ |24|), so the soft Wilson lines can be suppressed. This is why one only 
finds SlSn for this operator in the literature. For typical regulator choices the other gluons are simply 
not required to reproduce the infrared structure of the full theory result. 
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Here /„. ~ 131 MeV and the QCD field ip denotes the isospin doublet {if>M,il)W}. The 
coordinates satisfy (y — x) 2 = 0, which ensures that the path from y^ to x^ is along the 
light-cone, and Y(y,x) is a Wilson line along this path. In SCET we require the matrix 
elements of highly energetic pions, which therefore have collinear constituents. Boosting the 
matrix element in Eq. ( |32"D and letting y^ = yn^, x^ = xn^ we have 

(7r^ p \Cn,yKW(y,x)^ x \0) = -iUn- P 5 ah [ dz e «H**+(i-*)*] ^ z ), (33) 

Jo 

where = ^7sT fc / 'a/2, and £ n)2 is a collinear field with position space label zn^ . For our 



purposes it is more useful to use the operator with momentum space labels [13 



(K, P \Zn, Pl Wr b J(uj-P + )W^ n , P2 \0) = J ^e-^« p |^r^(y,-y)a,-,|0> (34) 



= —if 1T n-p5 a I dz5[co — (2z—l)n-p](j) 7! -(fjL,z), 
Jo 

where V± =V^±V. In Eq. (0) the delta function fixes u to the sum of labels picked out 
by the V+ operator. The combination picked out by V- is equivalent to (— P) acting on the 
entire operator, and using Eq. @ is fixed to the n-p momentum of the pion state. 

In some situations it is convenient to have delta functions which fix the labels of both 
W^£ n ,p and £, n , p W. In this case a useful field is 

= [8{o> - P) WjeS] ■ (35) 

Here i is the flavor index and will be omitted if the flavor doublet field is implied. Note 
that unlike the p in the label u on \ n w is not summed over. A matrix element with \n,uj 
fields is related to a matrix element like the one in Eq. (|34]) through 

(M n , P \Xn,^Xn^\M nip/ ) = 2«J(w_-fl-p_) (M n>p \ £ n , pi WT5 (UJ+ -P+)W^ n7 p 2 \ M n y ) , 

(36) 

where u± — to ± u' and p_ = p — p' . Thus, with the x notation the momentum conserv- 
ing delta functions become explicit. The factor of two appears from treating the u's as 
continuous variables, and in the final results cancels with a factor of 1/2 from a Jacobian. 

For inclusive processes such as DIS it is the proton parton distribution functions for 
quarks of flavor i, fi/ p (z), antiquarks fi/ p (z), and gluons, f g / p (z) which are needed. The 
standard coordinate space definitions [^] are (y^ = yn^ 



fi/pV) = l^e-^y(p\^(y)y(y,-y)^ { H-y)\p) _ _ , (37) 



spin avg. 



f g / P (z) = -J- /|fe- 2 -^nX<p|^ A (y)F afe (y,-y)Gt(-y)|p> 



spin avg. 



and fi/ p (z) = -fi/ p {-z). Here G° x (y) is the gluon field strength, Y(y, -y) and Y ab (y, -y) 
are path-ordered Wilson lines in the fundamental and adjoint representations, and \p) is the 
proton state with momentum p. In SCET these distribution functions can be defined by the 
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matrix elements of collinear fields with collinear proton states 



\Y.(Mx$u$xZ'\Pn) = 4n- P JdzS(uJ)S(u + -2zn-p)f i/p (z) (38) 

spin " 



— An-p dz 5(u)-)5(uj + + 2zn-p) fi/ p (z) , 
Jo 

J2(Pn\^ [B^B n /]\ Pn ) = -^p. £dz5(u„)5{u + -2zn-p)f g/p (z), 



spin 



where u± = uj ± uj' , and Bfc^ = n u {Q nuJ ) v ^ with the collinear gauge invariant field strength 



i9n*T = ~ [S(tv ~ P)W^VZ + gA^ q , W: + gAl q ,\w\ . (39) 



Both operators in Eq. ( p8|) are order A 2 since £„ iaJ ~ B^ u ~ A. Note that the matrix element 

of a single operator {Xnt^Xnt') contains both the quark and antiquark distributions. This 
is due to Eq. (|), from which we see that for uj = uj' > {uj = uj' < 0) this operator reduces 
to the number operator for collinear quarks (antiquarks) with momentum uj. 

Processes other than DIS sometimes depend on more complicated distribution functions. 
In deeply virtual Compton scattering (DVCS) we will need to parameterize the matrix 
element of an operator between proton states with different momenta. In terms of QCD 
fields ip the nonforward parton distribution function (NFPDF) defined by Radyushkin in 
Eq.(4.1) of Ref. [27] are (up to a trivial translation) 

(p 1 , a'\i>®(y)Y(y, -y)^{-y) \p, a) (40) 



+ h(a', a) [ dz [ e in ^ 2z -Ov xf {z; t) - e -™-p{2z-Qy xf ( Z] t )] 



where t = (p — p') 2 , and ( = 1 — n-p'/n-p. Here e(er, a') and h(a,a') are matrix elements 
which respectively preserve or flip the proton spin. They are defined in terms of the proton 
spinors 

e(a', a) = u(p', a') ft u(p, a) , h{a', a) = u(p', a') u(p, a) , (41) 

where m p is the proton mass. The NFPDF for gluons is similarly given by [^7| 

nX<P>1 G? (y)Y ab (y, -y)Gl(-y)\p,a) (42) 
= ^ e{a', a) [ dz [^ 2z ~Oy + e -iMa*-Ov] q 

+ HLl h ( a > t a ) f dz [ e «-P(2*-C)i/ + e -ifi-P(2*-C)tf] jq^ ; t ) . 

2 /n 
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In SCET the definition of the NFPDFs in terms of collinear fields are 

-i 



{plAx%^xZ\Pn,°) = 25{u-+n-pQ dz (43) 

Jo 

x {e(V, or) [6(u+-(2z-C) n-p) ^\z; t) - 8(u++(2z-() n-p) ff(z; t) 
+h{a', a) [5(u+ - (2z~C) n-p) K.f (z; t) - 5(u + + (2z-() n-p) Tcf (z; t)] 

(p' n , a'\Tr [B^B n /] \p n , a) = 6(u_ + n-pQ j\z 

x je(er', cr) 5(u + -(2z-C)n-p)F%(z;t) + h(a',a)5(u + -(2z-()n-p))C 9 .(z,t)\ , 

where the spinors in e(a', a) and h(a', a) are two component effective theory spinors, so that 
u — > u n where 7jiu n = 0. Note that for p' —> p both £ — > and t — > 0. In this limit the 
NFPDFs reduce to the standard PDFs: 



lim ^ = fi/ p (z), lim Tl = fi/ P (z) , \im F[ 9) = zf g/p (z). (44) 
p— tp ' p—tp' p—tp' 



D. Symmetries for collinear fields 

In this section we discuss spin and discrete symmetry constraints on operators involving 
collinear fields. 

The possible spin structures of currents with fields is restricted by the fact that they 
have only two components, ^ n>p = 0. The four most general spin structures for currents 
with two collinear particles moving in the same or opposite directions are 

ln,p' Tl in,p Ti = {ft , $75 , $7±} , 

£n,p> T 2 £n,p T 2 = { 1 , 75 , 7j£ } • ( 45 ) 

Other choices for Fx and T 2 either vanish between the fields or are related to those in Eq. (|45*|). 
This result can be expressed in a compact way by the trace formulae 



W r ^ p , r = |ir m - ^Tr[^ 75 r] - ^Tr[^r] , 

£ n j r Up , T = ^Tr[#r] + ^Tr [ 75 #r] + ^Tr[ 7/ ^#r] , (46) 



which reduce a general T to a linear combination of the terms in Eq. fl45|). For instance, 
it implies that 2i^ n a^ u i n = n v ^ n ^y^ n - n^ n fYxin, and £n7±75£n = i^inlv^n where 
e vv = e^ ual3 n a ni3/2. Furthermore, each of the two components of £ ra and also £ s can be 

chosen to be eigenstates of their helicity operators, h = p ■ S with eigenvalues ±1/2. For 
these fields h is equivalent to the chiral rotation, h = 7s/2. The structures in Eq. ( f£5| ) split 
into two classes depending on whether they conserve or flip the helicity 

chiral even: £,n, P >{$, $75}£n,p , ln, p Hi£n,p , (47) 

chiral odd: £n, P '$l±£n, P , £n, P '{l, 75 Rn,p ■ 
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Since gluon interactions in QCD preserve helicity, integrating out hard QCD fluctuations 
results in effective theory operators with the same helicity structure as the original operators 
at leading order in A. 

The presence of labels on the effective theory fields makes their transformation properties 
under the discrete symmetries C, P, and T slightly different than in QCD. For example under 
P or T we have n^nso these transformations relate collinear fields for different directions. 
Under charge conjugation, parity, and time-reversal the collinear fields transform as 

C~Hn, P {x)C = -[a,- P (*)Cf, C-^ iP (x)C = -[A^ p (x)f, 

P- 1 ^ p {x)P = 70 £*${xp) , P- l A^ p {x)P = g^KsM , (48) 

T- 1 e„, p (x)T = n nS {x T ) , T- l A» iP {x)T = g, u Al P {x T ) , 

where C~ lr y^C = — 7J and T = 7 5 C, while if x^ = (x + ,x~,x ± ) and p^ = (p + ,p~,p ± ) then 
p^ = (p~ ,p + , —p±), x p = (x~ , x + , — ar 1 ), and Xj, = (—x~,—x + ,x ± ). The transformation 
properties of W n can be worked out using Eq. ([KJ), for instance C~ l W n C = \W%\ T . 

The collinear effective Lagrangian fllOD is invariant under the transformations in Eq. (|48[) 
(adding the n «-> n terms). These symmetries also constrain the form of non-perturbative 
matrix elements. As an example, for a meson which is an eigenstate of C one finds 

(M n \^ p W n ^ 5 8(u - P + )W^ p/ \0) 

= (-l) c (M n \(CW^ p ) T ^5(uj - V + )(U P W n C) T \0) 

= (-lf(M^n, P 'W n fL l5 5(uj + V + )W^ n , p \0) . (49) 

For the isotriplet pion state (— l) c = +1 so combining Eq. ( f49"|) with Eq. ([34] ) gives 

(KjU Pl WT b J(uj-P + )W^ n , P2 \0) 

= —if^ ah n-p I dx 5[—uj — (2x — l)n-p](f) 7T (x) 
Jo 

= -ifJ ab n-p I dxS[u-(2x-l)n-p]<f> ie {l-x). (50) 
Jo 

Together with Eq. (0), charge conjugation therefore implies that <p w (l — x) — (j) n (x). 

III. EXCLUSIVE PROCESSES 
A. 7T-7 Form Factor 

The pion-photon form factor F nl (Q 2 ) is perhaps the simplest setting for factorization 
since there is only one hadron in the external state. The form factor is measurable in single- 
tagged two photon e~e~ — > e~e~n reactions. This process involves the scattering of a highly 
virtual photon and a quark-anti-quark constituent pair off an on-shell photon. The photon 
scatters the quark pair away from the incoming photon into a pion, so that 7*7 — > 7r°. The 
matrix element for this transition defines the 7T-7 form factor 

(7^)1^(0)17(^,6)) = iee v JdHe-^- z {^)\TJ^)J v {z)\Q) 

= -ieF^{Q 2 )^v P <rP u yq°- (51) 
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FIG. 1: Tree level matching onto Oj in the Breit Frame. The graphs on the left include u and d 
quarks. 

Here J M = ipQ^^ip is the full theory electromagnetic current with isodoublet field ip and 
charge matrix Q = r 3 /2+l/6, and — q 2 = Q 2 ^> Aq CD where q = p^—p^ is the virtual photon 
momentum. It has been shown that the form factor can be written as a one-dimensional 
convolution of a hard coefficient with the light-cone pion wavefunction ]28]. Here we show 
how this factorization takes place in the SCET. 

In the Breit frame q^ = Q(n tl —n tl )/2, the real photon's momentum is p^ = En^ ~ Qn fl /2, 
and the pion is made up of collinear particles with momenta n-pi — Q. The particles 
exchanged between the two currents in Eq. (Ml) have hard momenta and can be integrated 



out. At leading order in A the time ordered product of the two currents in Eq. fl5T|) matches 
onto a single operator in the effective theory. For simplicity we restrict ourselves to the 
tensor and spin structures that are relevant when the meson is a pion 4 



O 



Q ^ 



(52) 



0\ is of dimension 



where £„ iP is an isodoublet collinear quark field, and 2e^ v 

two, just like the time-ordered product in Eq. ([51]), and a power of 1/Q is included to 
make C(jU, V, V ) dimensionless. The time-ordered product in Eq. ([y]) is even under charge 
conjugation, so the operators in Eq. fl5"2"D must also be even. This implies CV 7 (/i, P, V ) = 



C nl (fi, —V , —V). The location of the Ws in Eq. (|52|) is fixed by gauge invariance, and V 
contains the spin and flavor structure 



(^(Sv^Q 2 ). 



(53) 



Aq CD we can also 



Since the offshellness of the collinear particles in the pion is p 2 
integrate out offshell modes with p 2 ~ QAqcd which come from soft-collinear interactions. 
For the collinear operators Oj, Eq. ( |PTD implies that factors of the soft Wilson line S n are 
induced. However, the location is such that SlS n = 1, so no coupling to soft gluons occur 
at leading order. The coupling of the collinear fields to usoft gluons can be simplified with 
the field redefinitions in Eq. (0). As discussed in section 11 A this moves all couplings into 
Oj, and using Y%Y n = 1 gives 



(54) 



Thus, usoft gluons also decouple. 



Note that a pure glue operator would not have the same isospin as the pion state. 
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In the Breit frame the pion momentum satisfies p% = E^n^ + 0(A), and comparing 
Eq. ( [5TD with the SCET matrix element i{7i^ p JO nl \0) , gives 

Y FAQ 2 ) = ^«IS^°)rc'(p,p^)^ )t^|o>. (55) 

Defining P± = ± P, the operator V- is related to P acting from the outside on the 
fields. Using Eq. (|6|) it can therefore be set equal to the momentum label of the state, 
V- = n ■ p n = Q. Suppressing this dependence we write C(P,P , //) = Ci(P + ,jj) leaving 

= ^ jduC^^^^lW^TS^ (56) 



Using Eq. fl34|) the remaining matrix identity in Eq. (|56|) can be written in terms of the 
light-cone pion wavefunction 

F ^(Q 2 ) = ^ JduJ dx5(u-(2x-l)2E v )C 1 (u, f i)(t> ir (x,fi) 

= tI / dxC^x-^Q^)^^). (57) 
v Jo 

This is the final result and is valid to leading order in A and all orders in a s . From Eq. (|50|) 
charge conjugation implies that 4> n (x) = 0^(1 — x) and Ciiu) = Ci(—u). Eq. (|57|) agrees 
with the Brodsky-Lepage |25j result that the form factor can be written as the convolution 
of a short distance function with the light-cone pion wavefunction. The SCET formalism 
gives a concise derivation of this result and defines the short distance function in terms of 
the Wilson coefficient of an effective theory operator. 

As an illustrative example consider the tree level matching onto C illustrated in Fig. [I]. 
Since the location of the PU's in O are fixed by gauge invariance, C([i,P,P ) can be de- 
termined by matching with W — 1. Expanding the full theory graphs to leading order 
gives 

i (Fig.l) = % W/^nV (t 75 ) (Q 2 ) (J--J-Y (58) 
2 \2 / \n ■ p n ■ p' J 

where we have dropped isosinglet terms, contributions with opposite parity to the pion, as 
well as those proportional to ^y 5 since ^£„ iP = 0. Comparing Eq. fl55| ) to Eq. (^) gives 



Q ( l l 

so that 



C(*P,P') = ^=(^ " i) + , (59) 



CiOi, u = (2x- 1)Q) = ^= (i + ^4^) + 0(a.(Q)) . (60) 



This result is again in agreement with Ref. p5[ , and the order a s (Q) corrections to this 



Wilson coefficient can be read off from the results in Ref. [29, pQ]. An identical analysis 



applies for operators with different spin structures such as the ones contributing to 7*7 — > p . 
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FIG. 2: Tree level matching onto Co s- The QCD graphs on the left plus the analogous graphs 
with the current on the bottom quark line are matched onto the collinear operator on the right. 

B. The large Q 2 meson form factor 

Another example of an exclusive process which can be treated in the effective theory 
is the classic case of the electromagnetic pion form factor at large space-like momentum 
transfer. For generality we consider in this section the electromagnetic form factors for 
arbitrary mesons (pseudoscalar P or vector V), defined as 

(P'(p , )\j»\P(p)) = F P (Q 2 )(p»+p' lx ), (61) 
(V(p',e')\j»\P(p)) = G{Q 2 )ie^ a ppy a e'P , 

(V'(p',e')\j,\V(p,e)) = F x {Q 2 ){e'* • e){p +p% + F 2 (Q 2 )[(e'* ■ p)e„ + (e ■ p')e';\ , 

where q 2 = —Q 2 ,q = p — p'. For simplicity we suppress the dependence of the form factors 
on the isospin of the two mesons. We will restrict ourselves in the following to the case of 
hadrons made up only of u, d quarks. The electromagnetic current is defined as usual by 
J[i — Q Qlfi Qj with charge matrix Q = diag(2/3, —1/3), which can be written in terms of 
the up and down quark charges as Q = (Q u — QdjTs/^ + (Q u + Qd)l/2. 

We will be interested in the asymptotic form of the form factors in the region with 
Q 2 ^> rrip, where it can be expanded in a power series in 1/Q 2 [P5[| . It is convenient 
to work in the Breit frame, where the momentum transfer has the light-cone components 
q = (q + ,q~,q±) = (Q, —Q,0). In this frame the meson momenta are p = (Q,mp/Q,Q), 
p' = (m 2 p /Q,Q,0), so the partons in the incoming/outgoing meson are collinear along the 
n^/n^ direction. 

The electromagnetic current in Eq. ( |BTD is matched in the effective theory onto the most 
general combination of operators constructed from collinear fields which are compatible with 
collinear gauge invariance. Operators such as the dimension three current 

[£nW n ]rc(p,r\P)\wk»], ( 62 ) 

can contribute, but only overlap with the asymmetric meson states with one energetic 
collinear quark and one usoft or soft quark. Often this overlap is referred to as the tail 
of the wavefunction contribution or the Feynman mechanism of generating the form fac- 
tor |TT], |52]| . There are other operators with significant overlap with more symmetric meson 
states (where all the constituents are allowed to be energetic). The leading such operators 
have the form 5 

[U Pl W n TWk %P2 ]C{fx,P,v\v,V^^ , (63) 



5 There are also gluon operators that can contribute when one or more of the mesons is a neutral isosinglet, 
however for simplicity these are not discussed here. 
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with C a dimensionless Wilson coefficient. As usual, collinear gauge invariance is enforced 
by the location of the W's in Eqs. (|62|) and fl63|). There is some argument about the relative 
size of Eqs. ( |52]) and (|63| ) in the literature |32|]. Often it is argued that the tail of the 
wavefunction is suppressed by an extra Aq CD /Q 2 pif , in which case the operator in Eq. ( |63"D 
dominates by two powers of Q. An analysis of the tail of wavefunction contributions has not 
yet been performed in the effective theory framework. Therefore, we choose to ignore the 
operator in Eq. QB"2]), and below only analyze the operator in Eq. fl6"3"|). We emphasize that 
we do not claim to have shown that this is justified by the effective theory power counting. 

There are two different structures possible for the operator in Eq. (^), and we write the 
general matching for the electromagnetic current as 



C (n, Uj) Jo(uj, //) + C 8 (/i, Uj) J£((i, Uj) , (64) 
where j = 1,2,3,4. The SCET currents are dimension-6 operators 

Jq = Xn^Xn,^ Xn^'Xn,^ ~ (T <-> V , UJl >2 «"> -^4,3) , (65) 
Js = Xn,u n TTa Xn,u 2 Xn,uX Ta Xn,u A - (IT <-> V , W 1)2 -^4,3) , 



where the x fields are defined in Eq. (p5|). In terms of the charge matrix Q, the spin and 
flavor structure is 

r<8>r' = {n v +n v )C^Q® 1 a 1 l\ . (66) 

The Wilson coefficients Cq^ can be computed in a power series in a s (Q). They are functions 
of n, Q, and the Uj which are the sum of momentum labels for gauge invariant products of 
collinear fields in the SCET currents. 

The currents operators in Eq. fl65|) are the most general allowed operators which are 
gauge invariant, transform the same way as J M under charge conjugation and satisfy current 
and helicity conservation. To see how these properties constrain the form of the allowed 
operators, we begin by noting that Eq. ( f45|) implies that T,T' = {1,75,7^} are the most 
general allowed spin structures. For massless quarks the electromagnetic and QCD couplings 
preserve helicity, whereas £n{l,75}£n cause the helicity to flip. Thus, only the structure 
Cn7x^n is allowed. Current conservation q v ' J v = 0, together with q u = Q(n u — n u )/2 implies 
J v oc (n u +n u ). Under charge conjugation J M — > — J M so the same must be true for the SCET 
currents. In the current operators, charge conjugation switches Ui <-> — cj 4 , co> 2 —^3, and 
T <-> T', as can be seen from Eq. (4~3). Thus, the second term in J~ U 8 is required to make these 
operators odd under charge conjugation. The operators J Q V % and the full electromagnetic 
current are invariant under a combined PT transformation. This requires that the Wilson 
coefficients are real. 

The operators j7o,s are responsible for the P„ —>■ P n transition, while the reverse tran- 
sition P n — ► P n is described by similar operators with n <-> n. Parity invariance re- 
quires the Wilson coefficients of these operators to be identical to Co$(u)i). Demanding 
hermiticity of the electromagnetic current in the effective theory then gives the relation 
C 0)8 (ui, u 2 , u 3 , U4) = Cq^u^UxiOj^ojs). Since the coefficients are real they must therefore 
satisfy Co l8 (wi, u 2 , u> 3 , w 4 ) = C 0)8 (a;2, u it w 4 , u 3 ). 

To compute the matrix elements in the effective theory, it is convenient to Fierz transform 
the four-quark operators in Eq. (|65]). This gives 



C J + C 8 J 8 = K + n v ) °jJj » ( 67 ) 

3=1 
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where 



3 3 [Xn^i-^jXn,^] [Xti^^-Xti,^] ■ 

The spin, flavor, and color structures are 

1 



(68) 



ri ® r; 



iQu - Qd)ie 3bc {r b <g> r c ) [ft ® i + $ 75 <g> ^ 75 ] 



(Q« + ® 1 + r a <g> r a ) + (Q u - Q d )(l ® r 3 + r 3 

while r 3j4 = T a r 12 and Tg 4 = T a r' 1 2 . The new Wilson coefficients are 



1) 



(69) 



Ci(/W) 



l{ 1 ~ W) c ^ Uj) + W Co{ ^ 

c c 



+ (^1,2 <-> -^4,3) , 
- (^1,2 <-> -^4,3) , 



(70) 



with similar relations for C3.4 which are also in terms of Cq^- 

A few general predictions follow from the form of the operators in Eq. (^). 6 For mesons 
with spin, only helicity conserving form factors appear, and furthermore no off-diagonal (e.g., 
P — > V) matrix elements are present at leading order in 1/Q 2 . These results agree with 
Ref. [p8| . We also see that the form factors between arbitrary meson states are determined 
at leading power by only two hard coefficients, C and C 8 . 

Now consider what factorization tells us about the matrix element of the operators in 
Eq. (p8|). For the decoupling of usoft and soft gluons we will follow section [LI A] . Integrating 
out offshell modes with p 2 ~ QAqcd induces soft Wilson lines S n and S n , while the field 
redefinitions in Eq. make all couplings to usoft gluons explicit in the operators. Together 
these give 



<jv _ r-(o) ytctr.c y Y (o) ] [WO) ytctp/c y v (°) ] 



(71) 



Consider first the color singlet currents j = 1, 2. Here the F's and 5"s all cancel using 
unitarity of the Wilson lines. Since the Af^ and A^ q gluons only interact with fields in 
the n and n directions respectively, collinear gluons are not exchanged between the n and 
n quark bilinears. Thus, the matrix element between states with particles moving in the n 
and n directions factors 



(n\Ji,2\n) = (n|xilr li2X wjo)(o|xilr; !2X i°L 2 |n). 



(72) 



Next consider the currents 3-$^ which have color structure T a £g> T a in Tj £g> T'j. In this 
case the usoft and soft gluons do not cancel, but can all be moved into one quark bilinear 
using the color identity Y^SU a S n Y n ® Y^StT a S n Y n = T a ® YislS n Y n T a Y^SlS n Y n . After 
this rearrangement it is clear that the (u)soft gluons and A£ and A^ q gluons only interact 
with the fields in one of the quark bilinears. Thus, the matrix element (n\Js^\n) factors, 



These predictions depend on the dominance of the operators in Eq. (p2|) over those in Eq. (|62[). 
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similar to Eq. ([721). For color singlet states, however, the matrix element of an octet operator 
vanishes identically since 



n 



y (0) tV 0) |o) = 0. 



(73) 



Thus, the effective theory currents do not contribute to the form factors at any order 
in perturbation theory 

Eq. (|72"|) shows that for arbitrary meson states factorization occurs. It remains to show 
that the matrix elements in Eq. ( [72]) are given by a two-dimensional convolution with the 
light-cone meson wavef unctions. To do this we consider the simple example of the 0~ + — > 0~ + 
form factor for the charged pion. It should be obvious that the same steps go through for 
other meson states. 

The symmetry of the pion wavefunction (j) n (x) under charge conjugation (x — > 1 — x) 
implies that only the J\ current contributes. Thus, 

FMQ 2 ) = ^ y*dw,- C7i0*, ^o^C?/) |5c^r lX JS, 4 |o><o|5^r' 1 xS*|^(p)> ■ (74) 



The required matrix elements can be obtained from Eqs. ([34]) and (|3~6D with {tt^) = ^(ti 1 ± 
in 2 ) I y/2. The momentum conserving delta functions fix uj\ — tu^ = n ■ p' = Q and cu 2 — uj 3 = 



n ■ p 



n ■ p = Q, while the uo = u\ + uj^ and cu' = to^ + cu 2 integrations can be done with the delta 
functions. This leaves 



F w± = ±(Q U -Q d )& JdxjdyTi 



x,y, fj)<j> w (x } ti)<h(y, n) 



where Ti(x,y) is defined in terms of Ci(ui, u 2 , 0Jz, u^) as 

Ti(x, y) = d (xQ, yQ, (y - 1)Q, (x - 1)Q) . 



(75) 



(76) 



The coefficients Cj(/j,, uij), and therefore also Tj(x, y), can be obtained at the scale /1 = Q 
by a matching calculation, as illustrated in Fig. 0. For this purpose, it is sufficient to compute 
the matrix element of the currents with free collinear quarks. To lowest order in a s (Q), only 
C$(u>j,{i = Q) is nonvanishing 



This implies 



C (uj j ,n = Q) = 



T 1 (x,y,fi = Q) 



Csiujj,^ = Q) = 4na s (Q) 



Q 1 



CU3 004, 



4vr« s (Q) 
9 



1 

xy (1 



1 



x)(l-y) 



(77) 



(78) 



Using the asymptotic light-cone pion wavefunction <j) n (x) 
with Ref. |25|, 



F n± (Q 2 ) = ± 



8nf 2 a s (Q) 
9Q 2 



-, 2 



x 



: 6x(l — x) we find agreement 
8nf 2 a s (Q) 



Q 2 



(79) 



The order a 2 (Q) corrections to Eq. ([77]) can be found in Refs. []33], |34], [35| 
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IV. INCLUSIVE PROCESSES 



A. Deep inelastic scattering 

DIS is a process which is both simple and rich in physics. As such it provides an ideal 
introduction to inclusive factorization in QCD, which we study from an effective field theory 
point of view in this section. The aim is to prove that to all orders in a s and leading order in 
A the DIS forward scattering amplitude can be written as an integral over hard coefficients 
times the parton distribution functions. This is done by matching onto local operators in 
SCET. The properties of SCET are used to show that matrix elements of the leading local 
operator can be written as a convolution of a hard coefficient with the parton distribution 
functions for the proton. 

The first step is to understand the kinematics of the process. The hard scale Q 2 = —q 2 
is set by the invariant mass of the photon, and x = Q 2 /(2p-q) is the Bjorken variable. 
In the Breit frame the momentum of the virtual photon is g M = Q{n ti — n M )/2, and the 
incoming proton momentum is p^ = n^n-p/2 + n^m 2 /(2n-p) ~ n fl Q/(2x) + n tl xm 2 /(2Q) 
up to terms ~ m 2 /Q 2 , where m p is the proton mass. By momentum conservation the final 
state momentum is P£ = q^ + p^, which gives an invariant mass P\ = (Q 2 /x)(l — x) +m 2 . 
Values f — x ~ Aq CD /Q correspond to the endpoint region where the particles in X are 
collimated into a jet, while values 1 — x ~ Aq CD /Q 2 correspond to the resonance region. 
We will consider the standard OPE region where 1 — x ^> Aqcd/Q so that the final state 
has virtuality of order Q 2 and can be integrated out. In contrast, although the incoming 
proton has a large momentum component in the n M direction it has a small invariant mass 
p 2 = m 2 ~ Aq CD , and therefore is described by collinear fields in the effective theory. 

Consider the spin-averaged cross section for DIS which can be written as 

rfff =2^^ L ^^ M >- < 80) 

where k and k! are the incoming and outgoing lepton momenta with q = k! — k, L^ v is the 
lepton tensor, and s = (p + k) 2 . The hadronic tensor W^ u can be related to the imaginary 
part of the DIS forward scattering amplitude: 

W^(p,q) = -lmT^(p,q), (81) 

7T 

T,u( P ,q) = \Y.(v\ f M\p), f lu ,(q)=i J d 4 ze^ z T[.Uz)J u (0)}, 

spin 

where for an electromagnetic current J M we can write 

T„,(P, q) = ( - 9** + Qj f-)Ti(x, Q 2 ) + (p, + (p v + |0 T 2 (x, Q 2 ) . (82) 

As explained above, the intermediate hadronic state has invariant mass P\ ~ Q 2 . There- 
fore, one can perform an OPE and match T^^(g) onto operators in SCET. All fields in the 
resulting operators are evaluated at the same residual space time point, however, the pres- 
ence of Wilson lines and label momenta make the operators nonlocal along a particular light 
cone direction. These nonlocal operators sum the infinite set of purely local operators of a 
given twist, however this is built into the formalism automatically. To match we write down 
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the most general leading operator in SCET which contains collinear fields moving in the 
direction, and enforce the condition from current conservation q^T^ = 0. This leads to 

t-f( ? Of + f) + ^±^(EO.'' + f), (S3, 

where 

of = &w£c®{p + ,P-)wi(®, 

O] = n^n v tx[W\G n y x W Cl{P+,V-)W\G n )ZW], (84) 

where i labels the flavor of the fermions and ig G% x = [W£+gA% q , iV x + gA x q ,]. The Wilson 
coefficients are dimensionless functions of V+, Q, and fi. As in previous sections we 
can separate the hard coefficients from the long distance operators by introducing trivial 
convolutions. This gives 

Of = fdu x du 2 Cf{u^uJ)[^^ X % 2 ], 

Of = - jd^d^q{u+,uiJ)\x[B!^BF*] , (85) 

where u)± = uj\ ± u> 2 , and B^ w = n v {Qn,u) u ^ with {Q n ^)^ x defined in Eq. ([39|). Next we 
factor the coupling of usoft gluons from the collinear fields using the field redefinitions in 
Eq. ([13]). The operator of has the structure in Eq. ( pT|) so the Y's cancel trivially, while 
for Of we find 

BZ = Y n B^Y^ (86) 

and the factors of Y cancel in the trace. It is easy to see that soft gluons also decouple 
using Eq. (|17D or by noting that there is no non-trivial soft gauge invariant way of adding 
soft Wilson lines S n to Of or Of. Under charge conjugation the full theory electromagnetic 
current J M — > — J M and therefore the operator T Mi/ — > T^. This implies relations for the 

effective theory Wilson coefficients since the operators Of must also respect this symmetry. 
Thus charge conjugation gives 

J dux du 2 Cf (w+, u_) xSU^xS^ = ~ jdui du 2 Cf (uj + , w_) Xn-^X^-^ 

= Jd^duj^-Cfi-uj^uj^xtML- (87) 
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In the second line we changed variable u\ — > —uj2 and uj 2 — > —oji which takes uj + - 
and — > Thus, to all orders in perturbation theory C^\— uj + , uoJ) = —Cj(u + ,uJ). 
This relates the Wilson coefficients for quarks and anti-quarks. Note that the above results 
are all independent of the collinear hadron on which DIS is performed. 

Next we take the matrix element between proton states. Using the definitions of the 



nonperturbative matrix elements given in Section [II Q , and picking out the coefficients of 
the tensor structures we find that the delta functions in Eq. (|38f) set u + = ±2Q^/x and 
uj- = 0. Since charge conjugation relates negative and positive values of uo + only coefficients, 
Cj(u) + , 0), with positive u + are needed in the formulae for DIS. Therefore we define 



Cj{2Qz,Q,Q,ii) 



(88) 



where here we have made the dependence on Q and /i explicit. Combining this with Eqs. ( jS2[) 
(0), (H), and (|38|) , gives the final result 



T 2 (x,Q 2 ) 



* 1 Ff) (I) + + h m (!) 



4x 

Q 2 " 



+ 



4#. 



x 
(Off 



2 7/| 



(off 



2x L 



i // '(^l-//?^ 



.r 



fo(0 + fi/ P (0] 



(89) 



where a sum over % is implicit. The hadronic tensor components Wi^{x,Q 2 ) = 
Im T lj2 (x, Q 2 )/tt and therefore are determined by the imaginary part of the Wilson coeffi- 
cients. The Wilson coefficients are dimensionless and therefore can only have a s (Q) ln(///Q) 
dependence on Q. This reproduces the Bjorken scaling of the structure functions. 

Finally, consider the tree level matching onto the Wilson coefficients shown in Fig. |3|. 
From these graphs only the quark coefficient functions C- can be non-zero and we find 



Q 2 ir6(z 



ImH, 







(90) 



where Qi is the charge of parton i. The vanishing of Imif^ a ^ tree level reproduces the 
Callan-Gross relation W x /W 2 = Q 2 /(4a; 2 ). 



B. Drell-Yan, pp l+l~X 

Next we will extend the DIS analysis to the Drell-Yan (DY) process: pp —>■ £ + £~X. 
Specifically we consider the Q 2 distribution, where Q 2 is the invariant mass of the lepton 
pair. Drell-Yan is more complicated than DIS because one has two hadrons in the initial 
state. In the center-of-mass frame the incoming proton and anti-proton move in opposite 
lightlike directions, and to prove factorization we use the fact that collinear modes in different 
lightlike directions can only couple to each other in external operators in SCET. We take 
the incoming proton to move in the n M direction and the incoming antiproton to move in 
the n M direction. The hard scales in DY are Q 2 and the invariant mass of the colliding pp 
pair s = (p + p) 2 . The lepton pair has an invariant mass Q 2 , and the invariant mass of the 
final hadronic state is 

^H-^-i)- (91) 
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where 



Q 2 Q 2 Q 2 (Q0 , 

t= — , xi = - , x 2 = — — . 92 

s 2p ■ q 2p ■ q 

We are interested in the kinematic region where p\ ~ Q 2 ', which implies that both x\ and 
22 are far away from one. As r approaches one the invariant mass becomes too small for 
the treatment given here to apply. However, the effective theory can be used to deal with 
this region as well. It is also possible to study the q± distribution, but this again requires a 
generalization of the discussion given below. 

The spin averaged cross section for Drell-Yan is 

3W d 3 h d*k 2 

dG ~ » uW (27r)3(2fc?) (2tt)3(2A; 2 ) ' [J6) 



where 



w " u = ^EE( 27r ) V4) (p+^-^-^)^i JM ( )i x )( x i r (°)i^) 

spins X 

= \J2 J d 4 xe- iq - x (pp\J»(x) J"(0)\pp). (94) 



spins 



The sum over spins refers to the initial hadron spins (the sum over final hadron spins is 
included in the sum over X). Integrating Eq. ( P3| ) over the emission angles of the final 
leptons one obtains 

da 2a 2 



iQ > 3Q,. j EM*!**. < 95 > 

" " spins 

where we have neglected the lepton masses and defined the operator 

W(t, Q 2 ) = -J^L 9(q°)6(q 2 -Q 2 ) jd 4 x e-<** J»(x) J M (0) . (96) 

As we discussed above in the region of phase space under consideration p 2 x ~ Q 2 , so these 
hard fluctuations can be integrated out. Operationally this means we match W onto local 
operators in the effective theory. We would like to show that the minimal set of order A 4 
operators that contribute to Drell-Yan are 

w - i fa c„K Q) [xi% $ -xSU] [xil i xfij 

+ ^^(^^^Tr^^lTr^B^] , (97) 



25 



FIG. 4: Tree level matching onto the operators in Drell-Yan. 

where the powers of Q are included to make the coefficients dimensionless. The operators 
displayed in Eq. ([57]) are just products of the operators that occured in DIS, so for these 
terms the decoupling of soft and usoft gluons occurs in a straightforward manner. To show 
that the operators in Eq. (R7J) are the most general set needed we must show that all other 
operators that are order A'Hjither reduce to these or vanish between the matrix elements in 
Eq. d95|) . For instance, A 4 operators also exist where field is contracted with a B^ w , 

field, or the color structures of the operators in Eq. ( j97j) could be arranged in a different 
way. 

We now give a general argument for why we can always rewrite an arbitrary operator in 
the form of Eq. ( |97D or show that it does not contribute to DY. All operators relevant for 
DY contain four order A collinear fields chosen from £ njP , B£ p , or B£ p . Furthemore, two 
must move in direction n and two in the direction n (other possibilities end up vanishing by 
baryon number conservation or because they involve a set of fields between physical states 
that can not possibly form a color singlet operator). For operators with 4 quark fields, 
Fierz transformations can always be made to arrange the fields such that those in the same 
direction sit in the same bilinear. Using as an example the operator with two collinear 
quarks in the n directions and two gluons in the n direction and leaving out the soft Wilson 
lines for the moment, the most general matrix element is 

(v v-\ y ' y m B-'^ B B,U Id v-) A ab ' AB (Qg) 

where a, b are quark colors, A, B are gluon colors, and a, ft are spinor indices for the quarks. 
^uv-ap i s some tensor that connects the indices in an arbitrary way. In the contraction of 
a, b and A, B there are two possible ways to make an overall color singlet, one where both 
the quarks and gluons are in a color singlet, and another where both the quarks and gluons 
are in a color octet. We will discuss both of these possibilities in turn. 

In the color singlet case, including the soft and ultrasoft Wilson lines is trivial, since 
using Eqs. (|i~7D, (|T^) , and (|A6|) we see that they cancel due to unitarity/orthogonality of the 
various Wilson lines in the fundamental/adjoint representations. Thus, there are no soft, 
usoft, or collinear interactions that connect the n and the n fields. As in previous sections 
this leads to a factorization of the matrix element in Eq. (|98|) , namely 

(Pnl xSlTx^f \Pn) (Pn\B^t B / \Pn) • (99) 

Since the proton spins are summed over, we can write (with the help of Eq. (]46|) ) 

(Pn| X { tTX [ Sf \Pn) OC 5 a \iU (PnlX^xtr W > ( 100 ) 

so that spin and color are summed over in the matrix element. Similarly the antiproton 
matrix element can be simplified to 

(Pn\ BZ A /BZ B : U \Pn) oc 5 AB gT(Pn\ Tr [B^B^] \p n ) . (101) 
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Here we used the fact that the matrix element is symmetric in /i and v, and that only the 
perpendicular index of the field is order A. Using Eqs. ( |100| ) and ( |101|) the original 
matrix element in Eq. fl98|) can be written as 

(PnPn\ ^n,W3^,U \PnPn) (102) 

oc Tr[A£ i\ ( Pn \ Tr [B$>B^] \p n ) (pj j& jfrg, \Pn) , 

where the trace of A is over spin and color, and just gives an overall constant. The final 
result in Eq. (|102|) is identical to the matrix element of the second operator in Eq. (|97j). 

If each of the n and n field bilinears involve color octet structures, then the soft and usoft 
Wilson lines don't cancel, since they don't commute with the SU(3) generators. However, 
one can use the color identity 

^n^liX SfiYn ® YftSftT S n Y n = T ® S n Y n Y^ S\T S n Y n Y^Sli (103) 

where each T x , S, and Y factor is in the appropriate representation of the color group. 
Eq. ( |103| ) moves all the soft and ultrasoft interactions between either the n or the n collinear 
fields. Thus, again the fields in one bilinear can not be contracted with fields in the other 
bilinear and the matrix element factors. However this time the factored matrix element 
vanishes. For the example discussed above, 

(Pn\X i Z a T C X i Z"\Pn}=^ (104) 

since an color octet operator vanishes between color singlet states. The same holds true for 
the matrix element of an octet gluon operator. 

An identical proof of decoupling goes through for the case of 4 quarks, where we again 
either have two color singlet or two color octet n and n bilinears. With 4 gluon fields we 
can either have the two n and two n fields coupled as singlets, or coupled in the same higher 
representation (an 8, {10, 10}, or 27). In the latter case the matrix element between color 
singlet states still vanishes so the proof for the 4 gluon operators also goes through in an 
identical way. 

Thus we have shown that the matrix element of an operator with an arbitrary contraction 
of indices either vanishes or can be written in terms of a product of a matrix element which 
is related to a proton pdf and a matrix element which is related to an antiproton pdf as in 
the example in Eq. (|102|) . This is the result we want. To see how the final formulae are 
derived note that we can write the matrix element of Eq. ( |9^) in the form of a convolution 

£ Y^{PnPn\ W \Pnpfi) = ^ f dUi C aA W i) (Pn\0*(uJ + , UjJ) \p n ) (p n \0^(u' + ,u'_)\p n ) , 
spins a,b 

(105) 

where u± = uu% ± u 2 and u>'± = ± U4. The operators here are the same as in DIS, with 
a = (i) for the quark operator, and a = g for the gluon operator 

Off(w+,a;_) 
Oi{u + ,uJ) 



1 r _ 



Q 2 



Q2 11 F«,^i fl f. J • 



(106) 
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Apart from the dependence on the labels, the Wilson coefficients in Eq. ( |105|) can also 
depend on the renormalization point [i and the kinematic variable Q. Using Eq. (R_8j) we see 
that the matrix elements in Eq. ( |106| ) set uo^ = oj'_ = and uj + = 2\fsz\, uj' + = 2^/sz2 where 
Z\ and Z2 are the convolution variables. Since all kinematic variables aside from Q 2 are 
integrated over in Eq. fl96|) the only other variable that can appear in the Wilson coefficient 
is the center of mass energy which produces the £ + £~ pair, namely 4s = uj + uo' + . Thus, the 
Wilson coefficients only depend on uo + uo' + = &sz\z 2 . Defining new coefficients 



H a \z lZ 2) = C ab (u; + u;' + = 4 S ziz 2 ,Q,n) 



(107) 



we can replace the matrix elements in Eq. ( |106| ) with parton distribution functions using 
Eq. (RSI) to obtain: 



T ^2(PnPn\W\p n p n ) 



spins 



— / dz\ dz 2 



■H { ^ {j \-z x z 2 ) fi/ p {z l )f j/p {z 2 ) + U/p^fj/pizz) 



+H m \ Zl Z 2 ) |/i/p(Zl)/j/p(z2) + fi/ P (Zl)fj/p(z 2 ) 

■ 22 H^Hz l Z2)f l/p {z l )f g/p {z2)--^H^' 

ZlZ 2 )f g/p {zi)fj/p{z 2 ) 



1 HrtHzi^L/JzAfi/Jz*) - A= 



2v^ 



2V? 



+ — H^( Zl Z2)f g/p ( Zl )f g/p (z2: 



4r 



-ZlZ2)fi/p{Zi)f g / p {z2) 
-ZlZ2)f g l v {zx)fjl p {z2) 

(108) 



This is the final convolution formulae for Drell-Yan and is valid to all orders in a s and 
leading order in the power expansion. At tree level the matching calculation shown in Fig. |] 
yields zero for all the Wilson coefficients except 



H^\- Zl Z2) 



2vrr 



Q 2 S(r - z x z 2 ) 



(109) 



where Qi is the charge of parton i. The coefficients H^' 9 (±ziz 2 ) and H 9 ^(±ziz 2 ) start at 
order a s (Q), while H^ 1 ' (z\z 2 ) and H"{z\z 2 ) start at order a 2 (Q). 



C. Deeply Virtual Compton Scattering, j*p —>■ "f^p' 

Next we examine deeply virtual Compton scattering (DVCS). To be more precise we 
examine the exclusive reaction j*p — > 'j^p, where the incoming photon is highly virtual, the 
final photon is either off-shell or real, and the incoming and outgoing protons have different 
momenta. The reason we have included this process in the inclusive section is that DVCS has 
the remarkable property that the nonperturbative physics is described by a so called non- 
forward parton distribution function (NFPDF). The NFPDF is a more general distribution 
function that reduces to the standard parton distribution functions (familiar from DIS) for 
some values of the momentum fraction, and behaves like a lightcone wavefunction (familiar 
from the pion examples) for other values. Deeply virtual Compton scattering was first 
studied in perturbative QCD in Refs |p_6|, KJ], p_8j, and proofs of factorization to all orders in 
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perturbation theory were later presented in Refs ^(J . In addition properties of NFPDFs 



were studied in Ref. [[27]]. Here we present a proof of factorization for DVCS based on SCET. 

As with the previous proofs it is important to understand the kinematics of the process. 
We take the incoming proton and photon momenta to be p and q respectively, with x = 
Q 2 /{2p- q) and q 2 = —Q 2 ^> Aq CD . The outgoing proton and photon momentum are p' and 
q' respectively, with > q' 2 > —Q 2 . It is convenient to define a parameter ( = 1 — n-p'/n-p, 
which measures the change to the proton's large momentum. Working in the Breit frame 
and neglecting contributions that are -C Aq CD /Q we have: 



Label Momenta Residual Momenta 
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+ <2Qm 2 p (n^+n^) 
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-1) n» + ^-n^-pf 





(110) 



Here the label momenta are order Q or QX, while the residual momenta are order QX 2 and 
depend on p 2 = m 2 and t = (p' — p) 2 = (p'l — ( 2 m 2 )/(l — (), which are both ~ Aq CD . The 
invariant mass of the intermediate hadronic state is (p + q) 2 ~ Q 2 {1 — x)/x just like DIS, so 
for 1 — x ^> Aqcd/<5 the intermediate state can be integrated out. 

We will proceed in a manner analogous to the analysis for DIS. The amplitude (up to an 
overall momentum conserving 5-function) is given by a time ordered product of currents: 

T^(p,q,q') = {p\v'\f^(q,q')\p,a) 

f^(q,q') = i J d*zj^*l 2 T\j^-zl2)J v {zl2)\ . (Ill) 

This time ordered product is contracted with a lepton tensor to obtain the amplitude. Now 
current conservation requires q^T^ = q' y T^ u = 0, however the DVCS T^ v is not symmetric 
under fi <-» v. For electromagnetic currents J M we have 

= -\ 9 ^ -wr 1 + & " ^ v v " T2 

+i/u T 3 + t»{ Vv - q -^f) T 4 + (p„ - T 5 + . . . , (112) 

where the functions Tj = Tj(x, (, Q 2 , t), and the vectors = q' tl —q fJ ,+p l j,(q 2 —q' 2 +t)/(2p-q) and 
= q^ — q^+Pfj,(q' 2 —q 2 +t) / (2p- q') are defined so that q-l — q'-i' — 0. In Eq. (|112| ) and below 

the ellipses denote spin dependent terms. For simplicity we will show how factorization is 

achieved for the spin independent contributions shown in Eq. (|112|) with the understanding 

that it is no more difficult to also include the other terms. 

It is convenient to define a parameter < a < 1, by q' 2 = — aQ 2 . The DIS hadronic 

time-ordered product is obtained in the limit p' — > p, where a — > 1 and ( — > 0. From 

Eq. ( |liq ) we see that 

,2. 



t Tfl 

( = x (l-a)+<D[— ,-f), (113) 
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so these parameters are not independent. Since the intermediate hadronic state has invariant 
mass 0{Q 2 ) we can match onto operators in SCET. Requiring q^T^ = and q w T ilv = 
for the order Q label momenta leads to 



T, 



9? 

Q 



w + 91 

1 Q 



+ - (n" + n") ( 



an" + n"l 



Of 

Q. 



+ 



(114) 



where the ellipses are spin dependent terms and the displayed operators are 



O 



tr[w\G n y x w q{v+,v-) w\G n )^w] . 

We have suppressed the dependence of the Wilson coefficients C(V+,V-) 



Of 



(115) 



on Q, a, and \i. 

The form of the operators in Eq. ( |115| ) looks the same as the DIS operators given in Eq. (^4|) . 
however the operators here are more general because the Wilson coefficients depend on a. 
In the limit a — > 1 the DVCS operators reduce to the DIS operators. However, since the 
field structure of the DVCS operators is identical to DIS several results follow immediately. 
For instance, the steps which factorize soft and usoft gluons and leave fields with superscript 
(0) are the same and are not repeated here, 



(<), 



(o)(0 



Of 



(0) 



(116) 



-tv[B^ q{v + ^) {B, 

The restrictions on the DVCS Wilson coefficients from charge conjugation are the same as in 
Eq. (|87|), Cj(V + , V-) = —Cj(—V + ,V~), however because p ^ p' this is not simply a relation 
between quark and anti-quark Wilson coefficients. The way in which DVCS is unique is that 
the matrix elements involve nucleon states with different momenta. This is what leads to 
results in terms of non-forward parton distribution functions. 

The definition of the NFPDFs are given in Eq. (^), and can be used along with the 
relations above to obtain expressions for the Tj in terms of the NFPDFs. Before we give 
this result we note that the Wilson coefficients depend on the operators V± which become 
the variables u± after introducing trivial convolutions and the Xt»m fields in Eq. (p5|). The 
delta functions in Eq. (13) then set oj- = —Q(/x and oj + = ±Q(2£ — ()/x, where £ is the 
convolution variable. Note that (/x = 1 — a, and just like DIS it is the combination £/x 
which appears. Since charge conjugation relates the Wilson coefficients for uj + > and 
uj + < it is convenient to define 

H j ({;/x)=C j (Q(2£/x-l+a),Q(a-l) : Q,a,Li), (117) 

where in the last three arguements we have made the dependence on Q, a, and /i explicit. 



Combining Eqs. 
Ti = 



e(a',a) 



, ( |116p , and ( |117| ) then gives 
i 



To 



2Q 

2 <l + a, 



<%<H\ 



(0 



X 



+ 



2x I 



2(1 + a) H% 



o 

i 

X 



2(1 + a)H 2 



X 



-Hi 



(0 



1 

2x 

i 

X 



H\ 



-H{ 

= 0, 



X/ J 
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which are the final convolution results valid to all orders in a s and leading order in the power 
expansion. The structure functions 23,4,5 vanish since the vectors £ M = = at leading 
order in the power expansion. The terms with ellipses are for the spin flip terms and involve 
the NFPDF /C defined in Eq. (|43|). The results for these terms have a similar form to those 
in Eq. flTTgp . 

Finally we match at tree level. The tree level diagram in QCD is the same as in Fig. |3| 
except the outgoing photon and proton have momenta q' and p' respectively. Only the quark 
Wilson coefficients are nonzero at tree level. We find 

which gives 

Since Ti,^ = at tree level, DVCS also obeys a Callan-Gross relation. 



-,Q,a) = e 2 Qf 
-,Q,a) = 0, 



2Q 



2Q 



2Q + uo + — oj_ 2Q — uo + — oj_ 



(119) 



\x 

) a 

X 



-e 2 Ql 



0. 



(120) 



V. CONCLUSION 



What we hope we have demonstrated here is the power of effective field techniques in 
the context of factorization for hard scattering processes. The explicit separation of modes 
and the implementation of gauge invariance for these modes greatly simplifies seemingly 
complex problems. What is normally accomplished by diagramatic Ward identities and 
induction techniques now falls out as a consequence of the gauge symmetry of operators in 
a low energy soft-collinear effective theory. 

As we have emphasized the factorization formulae derived in this paper are not new. The 
purpose here was simply to extend the formalism introduced in 0, |TUL [TT|, [TJ, |T^| to cases 
with back-to-back collinear particles, and apply these ideas to more general processes than 
previously considered. Furthermore, the factorization proofs presented are perhaps simpler 
than those previously given (certainly they are more concise). We believe that within the 
confines of the SCET more difficult, and unresolved problems can be addressed, such as 
power corrections in cases without an OPE, and proofs of factorization for more complex 
processes. 
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TvDe 


Momenta (+, — , _L) 


Fields 


Wilson lines 


onshell 


collinear-n 


rf~(A 2 ,l,A) 




w n 




collinear-n 


p£~(l,A 2 ,A) 




W n 




soft 


~ (A, A, A) 


9s, ^ 






usoft 


£^~(A 2 ,A 2 ,A 2 ) 


9ns, Aus 


Y Y- 


offshell 


p = pi + p 2 

p = Pl + q 
P = P2 + q 


(1,1, A) 
p»~ (A,1,A) 
p"~ (1,A,A) 


, A n M 


A>n : Afj 
" n ' n 

VF- X 



TABLE I: Summary of the onshell modes discussed in section HA, and the auxiliary fields intro- 



duced to represent the offshell fluctuations that are integrated out in this appendix. 



APPENDIX A: FACTORIZATION OF SOFT AND COLLINEAR n & n MODES 

This Appendix discusses the simultaneous factorization of the soft (A, A, A) modes, n- 
collinear (A 2 , 1, A) modes, and n-collinear (1, A 2 , A) modes. These three classes of modes 
can not interact with each other in a local manner and therefore do not couple through 
the SCET Lagrangian. However, they can couple in a gauge invariant way through external 
operators and currents. These interactions in currents are built up by integrating out offshell 
fluctuations with p 2 S> (QX) 2 . For the special case of factorization of soft from n-collinear 



modes this was shown in detail in the Appendix of Ref. [12]. There it was shown that 
integrating out certain modes with offshellness p 2 ~ Q 2 A causes the Wilson lines W n and S n 
to appear in operators in a gauge invariant way. Here we will extend this approach to the 
factorization of modes for cases involving two classes of collinear particles. For simplicity 
we restrict ourselves to the case where the original operators involve only collinear quark or 
gluon fields. This type of factorization was used for the pion form factor example discussed 
in Sec. III.B and the Drell-Yan process presented in Sec. IV.B. 

The basic idea is to first match onto a Lagrangian with couplings between onshell and 
offshell modes that give all order A diagrams. The offshell modes (with p 2 3> (QX) 2 ) can 
then be integrated out, so that all operators are expressed entirely in terms of the onshell 
degrees of freedom. In table | a summary is given of the three types of offshell momenta that 
are induced by adding soft, n-collinear, and n-collinear momenta. For each type auxiliary 
quark and gluon fields are defined, and for convenience momentum labels are suppressed 
in this section. For example, the ipQ quarks are created by the interaction of a n-collinear 
quark with an n-collinear gluon, whereas the quarks are created when a collinear quark 
£ n is knocked offshell by a soft gluon. For the field tJ)q we write ipQ = ipn + W , where 

ip® — z^'ipQ an d ipn = \TjitJiipQ. Then we have ifcipn = $W = an d = = 0- 
Various Wilson lines are also required and are listed in the table. It is convenient to define 
a generic Wilson line L[a, A] along direction a with field A by the solution of 

(a-V + ga-A)L[a, A] = 0. (Al) 

With this notation the on-shell Wilson lines are W n = L[n,A n ], W n = L[n,A n ], S n = 
L[n, A s ], and S n = L[n, A s ]. (Recall that the subscripts on W and S mean different things.) 



32 



The Wilson lines involving offshell fields that we will require are 

X n = L[n, A Q +A*+A n ] , X n = L[n, A Q +A*+A n ] , (A2) 
W* = L[n, A* + A n ] , W* = L[n, A* + A n } , 

S* = L[n, A*+A s ] , = L[n, A* + A s ] . 

Below we discuss the results which allow us to integrate out offshell fluctuations. The 
structure of the auxiliary Lagrangians and construction of their solutions are very similar 



to the case presented in Ref. ||12|| , to which we refer for a more detailed presentation. 

From table | we see that adding n and n-collinear momenta gives p 2 ~ Q 2 , whereas adding 
soft and collinear momenta gives p 2 ~ Q 2 X. Loops that are dominated by offshell momenta 
only modify Wilson coefficients and not the infrared physics. Therefore, to determine the 
structure of SCET fields in an operator it sufficient to integrate out the offshell fields at 
tree level. For convenience we can integrate out the fluctuations starting with those with 
the largest offshellness. Recall that we only wish to consider offshell propagators connected 
to external operators. A subtlety for quarks is that distinct auxiliary fields are needed for 
the incoming and outgoing offshell propagators. However, the solution for the outgoing field 
turns out to be the conjugate of the incoming field, so to avoid a proliferation of notation 
we simply denote the outgoing terms in the Lagrangian by +h.c, and present a solution for 
the incoming fields. Finally, note that for the gluon field Aq the fields A n , A n , A*, A* , and 
A s appear as background fields while for the fields A* and A* it is A n , A n , and A s that 
appear as background fields. 

The terms in the auxiliary Lagrangian involving the p 2 ~ Q 2 fields are 

= ^9niA Q + A^ + A n )^ + ^ n )+^[n-V + gn^A Q + A^ + A n )}^ 
+(n «-> n) + h.c. 

+i tr ( i W Q + 9 A Q > iD Q + ^q] 2 } + ^{[iVt,* A Q»] 2 } ■ (A3) 

where iDq = \n^\P + gn-(A^+A n )} + ^n^[P + gn ■ (A^+A n )}. The solution of the equations 
of motion for these modes are 

v# = (x n -i)(rt + q , $ = (x n -i)(^ + e ft ) , 

XlX n = W*wP. (A4) 

(In addition to the last equation a constraint on the components u-Aq and u-Aq also comes 
from the gauge fixing term, but will not be needed.) The terms in the auxiliary Lagrangian 



involving the p ~ Q X fields are [12 



£f ux = tP^gn-(A^ + A s )^ n + ^[n-V + gn-(A^+A s )]^ + (n^n)+h. C . (A5) 
+ ^r{[iD^ x +gA^,iD^ x +gA^) 2 } + ^r{[z^ x ,<J 2 } + (n <- n) , 

where iD^ x = \n^\P + gn-A n ] + ~n M [n-'P + gn-A s ]. The solutions for these modes are 

V£ = (#-l)£„, s^w* = w n si, 

It = (S*-m , S^W* = W n St . (AQ) 
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Together Eqs. and ( |A6|) can be used at leading order to eliminate the fields representing 
offshell fluctuations with p 2 ^> (QX) 2 . For collinear quarks this leads to the rules in Eq. (|TTD. 
Note that we did not need to use the gauge fixing term to resolve the ambiguity in the implicit 
solution for the fi-A and n-A auxiliary fields. 

As an illustration of these results we discuss the soft-collinear factorization for the pro- 
duction of a q n qn pair with a large invariant mass Q 2 . This process is mediated in the full 
theory by the electromagnetic current J = ipYip (Y a color singlet). This current will match 
onto a current in SCET that is built entirely out of onshell fields. Using the results in this 
appendix this current can be systematically derived. To start the quark field in J matches 
onto £ n plus all possible fields which the auxiliary Lagrangian can create starting from £ ra , 
so 

J^{L + ^ L n +^)T{u + ^k + ^)- (A7) 

Integrating out the p 2 ~ Q 2 fluctuations with Eq. ( |A4j ) and inserting a hard Wilson coeffi- 
cient C which depends on label operators turns Eq. ( |A7| ) into 

(£, + ft)Xl CY X n (C- n + V£) = {L + 4>nW* CY wFfa + ij>i) . (A8) 

To construct the first operator we used the equations of motion for ip® and ip®, and in the 
second operator we used the equation of motion identity for the gluons in X n and X n . In a 
similar fashion we can now integrate out the p 2 ~ Q 2 X fluctuations with Eq. (|A6|) to give 

LS^W* CY W^S^n = LW n Si CY S n Wkn ■ (A9) 

The operator on the right is the final result used in Eq. (|23J), and is soft, collinear, and usoft 
gauge invariant. It should be obvious from this example how the equations of motion in 
Eqs. ( |A4| ) and ( |A6| ) can be used to determine the factorized form of a general leading order 
operator. 
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